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Abstract. We study the diffusion of a particle on a random lattice 
with fluctuating local connectivity of average value q. This model is a 
basic description of relaxation processes in random media with geometrical 
defects. We analyze here the asymptotic behavior of the eigenvalue 
distribution for the Laplacian operator. We found that the localized states 
outside the mobility band and observed by Biroli and Monasson, in a 
previous numerical analysis are described by saddle point solutions that 
breaks the rotational symmetry of the main action in the real space. The 
density of states is characterized asymptotically by a series of peaks with 
periodicity 1/q. 
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Diffusion on random graphs can be a useful problem for studying 
relaxation processes in glassy systems in general. Usually, the disorder 
arises from a random potential, impurities, but a random geometry can also 
play this role One can visualize a diffusion of a particle on a random 
graph as the relaxation of a disordered system out of equilibrium on a 
complicate energy landscape. For example, this relaxation in random Ising 
magnets can be identified with diffusion on the vertices of a hypercube in the 
configuration space, and the edges correspond to the energy paths that connect 
one configuration to another j^. Numerical simulations jHl E] in Ising spin 
glasses of the order parameter q(t) = [< Si(t)Si{0) >]av, where < . . . > is the 
thermal average and [. . .]av the average over disorder, show that this quantity 
follows a Kohlrausch law similar to a "stretched" exponential exp(— (t/r)'^) 
with l/3</3<lin the region just above the spin glass phase. This kind of 
non-exponential relaxation is typical of many glassy systems |2l E] , unlike 
the usual exponential behavior with only one relaxation time. The coefficient 
f5 varies with temperature from 1 in the high temperature phase, to 1/3 at the 
glassy transition. In between, there seems to be a phase of localized states, or 
Griffiths phase, where f3 is in between from 1 to 1/3. The value 1/3 seems to be 
universal in many experimental systems and has been observed in spin glass 
Euo.4Sro.6S jZj, at the glass transition of the crystaUine fast ion conductor 
Na /3-alumina [H], or molten salts CaK(N03) [Oj. The question about the 
universality of the Kohlrausch law has been raised, and exact dynamics on 
ultra-metric spaces [10] show that it depends on the scaling of the barriers 
with the distance between different states. It could be algebraic (linear 
dependence with the distance) or Kohlrausch like (logarithm dependence). 
Random graph models can also be a useful tool to study propagation of sound 
waves in granular and random media fn\ I12j. since diffusion, quasi- diffusion 
and localization of sound waves are similar to diffusion on a random network, 
the edges being the connections between the particles making up the medium. 
Frequency response for different amplitudes shows a region of coherent 
propagation plus a quasi- diffusive regime. Numerical and simplified models 
of a granular medium for small vibration amplitudes show modes that are 
extended and localized in space, like a particle moving on a random medium. 

In order to simulate the glassy systems discussed above with short-range 
interactions and also percolation problems, Viana, Rodgers and Bray (VRB) 
studied a simple model of diffusion based on sparse random matrix ITH] . 
Their model has the property to have long-range interactions with dilute 
coordination numbers, so that in average the coordination number is finite 
and we may expect to have a closer view of some experimental material. The 
spectrum of eigenvalues (all positive) spreads over a continuous band, and for 
the particular case of the Cayley tree model, are bounded between mobility 
edges, with a gap from below (see also JE])- Recently, Biroli and Monasson 
(BM) py studied the same model numerically, and showed that localization 
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effects arise outside tlie extended region, witfi peaks at some regular intervals 
where the inverse partition number is high, showing that these peaks are 
resonance for localization. The existence of a Griffiths region was analyzed 
by Rodgers and Bray [Ej who found a tail distribution for large eigenvalues 
outside the mobility band (with a random matrix containing only elements 
-1, and 1) but there is no peak structure in their analysis contrary to BM's 
numerical result. Their saddle point solution is invariant by rotation in both 
replica and real spaces and they do not discuss the possibility of the rotational 
symmetry breaking in at least one of these spaces. We will analyze in this 
paper the possibility of a real space symmetry breaking for the site dependent 
fields, unlike the method of self-consistent field equations developped in [T7]. It 
is clear that in order to describe the localized regime, we have to take account 
of this kind of solutions. Our motivation is therefore to find the asymptotic 
solutions in the region of the density of states outside the mobility band for 
the VRB's model. The original model of relaxation in a geometrical disordered 
system consists of points connected randomly by bonds Mij = —1/q with 
a probability equal to q/N and Mij = else. This is an infinite range model, 
but the average coordination number q is finite and we may expect that the 
dilute exchange interactions make the model similar to a short range model 
|18j . We follow the references PP and PH] for the notations. On a random 
lattice, a particle performs a random walk from one point to another. Let 
Ci{t) be the probability for the particle to be on the site i at time t. Then the 
master equations describing the time evolution of these amplitudes are 

^ = -E^.-9 (1) 

where the probability for the symmetric elements Mij is 



for i j. The diagonal elements are equal to Ma = —J^j^iMij. This 
ensures that the eigenvalues of the matrix M are all positive since the quadratic 
form 



MijXiXj = -1/2 ^ Mij{xi - Xjf 



is always positive ^3]. Let Aj be an eigenvalue of M, with Vj the 



corresponding eigenvector, then the general solution of equation is 
c^{t) = E(V,),(V,)fcCfc(0)exp(-A,t). 
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If the particle is on the site m at t = 0, i.e. Cj(0) = 6i^rn, then the 
probabihty that after a time t the particle is found on the same site is equal 
to Cm{t), and we can define an average probability of return to the origin at 
time t after summing up over all the possible origin points: 



fit) 



j 



-X,t) 



P{\)exp{-\t)d\, 



(3) 



where [. . .]av is the average over the link configurations. -P(A) is the 
probability distribution function for the eigenvalues of the symmetrical matrix 
M. In figure 1 we have solved numerically P by diagonalizing 4500 different 
samples of random matrices M with the parameters = 800 and q = 20 (see 
also jl]). Then we made an histogram of all the different eigenvalues found 
and normalized the distribution. Following the results from the partition 
number of reference the distribution presents a mobility band between 
approximatively A = 0.5 and A = 1.7 plus localization edges outside this band 
and containing peaks. The distribution P{X) can be expressed as a functional 
of replicated fields 0f . Indeed, we can write 



P(A) 



tr5(A-M) 



-Im— trG'(A + ie) 

TlN 



where we introduce the Green function G{X + ie) = 1/(A + ze — M). More 
precisely: 



P(A) 



1 d 

-Im— — tr log( A + ie - M) 



— Im 



1 d 



log det(A + 26 - M) 



TiN dX 

The determinant can be rewritten as a Gaussian integral over fields (pi, 
1 N: 



P(A) 



2d r ^ ( i 



(4) 



Asymptotic behavior of the density of states on a random lattice 



5 



The average over the hnk configurations is performed by using rephcated 
fields 0f , and we obtain (see [Ej for the details): 



P(A) 



Im- 



2 d 



\ogZ 



lim Im —-^ 



lim Im \Z^' 

n^o irnN dX ^ 



exp 



ia2 



exp 



2 



1 - 



+ 



N N 



t>lf 



lim Im — 

n^o TrnA^ OA 



n#rexpS(A,{0r})- 



(5) 



Now we want to study the asymptotic behavior of the distribution P by 
evaluating the saddle points of S in the limit A ^ oo, while we keep large 
but fixed first. It seems very similar to studying the saddle point solutions 
of the extensive function S in the large N limit, but the difference is that we 
are studying the asymptotics of the distribution which does not describe the 
mobility part of the curve, so the corresponding saddle points are different. 
Moreover we are looking for solutions which can break the rotational space 
symmetry, since we are interested in the localized regime where strong or low 
connectivity may contribute to the distribution. A particle localized in one 
region would favor the amplitudes of the sites inside this region. The extrema 
of the functional S in equation are given by the set of equations 



(A + ^e)C 



-E- 

Nyi-q/N- 



exp 



g/iVexp(-^E^ 



.(6) 



The simplest solution is when all fields have the same value on every site 
0" = 0", which directly leads to 0" = 0" = 0. This gives a Dirac function 
centered at A = 1 and this is not an asymptotic solution for our problem. The 
next step is to take 0" = 0" everywhere except on one site iq where 0^^ ^ 0". 
We obtain 2 equations (up to the order 1/A^): 



(A + ie)<Pl 



(A + le) 



N -I 
N 



(j^l - 0") exp 




In the large N limit, only the field 0" does not vanish 




(7) 



-0"/(iV-l)), 



with the length of the vector 



" ) satisfying the equation: 
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0m = H (Pio = log (A + ie) + Anqm, and (p^^^^ 




where m is an integer. There is therefore an infinite set of non zero 
solutions in the complex plane. In the following we choose the site located at 
io = 1 for simplification, since the solution is invariant for the — 1 other 
sites. For each integer m, the number of all possible vectors with module 0^ 
is equal to N2t[^ /T{n/2) ~ Nn in the limit n — >■ 0. Given the solution (jH)). 
the value of the action is equal to 



S„, = -q{\ + ie) log(A + ie) + q{\ + ie) - q + 2i7rg(A + ie)m (9) 

then P falls exponentially with Alog(A), and the corrections to the 
exponential are given by computing the fiuctuations around the saddle points. 
Unfortunately the matrix of the second derivatives of S at the site = 1 has 
many zero eigenvalues and we need to consider the next order. Indeed, we 
have 



d^S _ X + ie o d^S i 



Tt^cP, i ^ j, 



90? 90? q 90f90f N 

___ = i(A + - 1)5^^, (10) 

The first matrix in (fTUI) has (n — 1) zero eigenvalues, and we therefore 
have to take account of the third derivatives at the site zq: 



X + ie 



')°5p^ + 01 (^07 + (pJSaf: 



b{d(Pl q 



+ (11) 

The fiuctuations around the saddle point solutions 0^ lead to the following 
integrals for main contribution at large N: 



te- LJi^x^x 



±_ sr _ A±i£^"^"^^^^ M ^(A + ^e) . ^ 



^x^x" - (piXiXiXi H —5 (f)i(f)i(f)iXiXiXi + 
2iV^ ^ 2q Qq^ 
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The Gaussian integration over the variables xf gives (27[ /i(\+ie))"'^^~^^^'^ . 
In the following we will set e = 0, since it is not essential for the remaining 
calculation. It is useful to work in the basis where the quadratic terms 
are diagonal. We first set 0° = 4>mf°' where {ip^ , . . . , (p"-) is a real unit 
vector, and (pm a square root of (jHl). The matrix ipa'Pis has one eigenvalue 
unity with eigenvector = {(p^ , . . . , ip"') and n — 1 eigenvalues zero with 
eigenvectors V^, P = 2, ..n satisfying J2a ^aVa = 0- We then use new variables 
Va = J2f3VpXi that diagonalize the quadratic terms. The matrix Vj^ satisfies 
*V = V~^, so that = J^f^VaVi^- After some algebra, can be finally 
written as 



A , ^-^ o iX 



-^<t>rr.yiJ2yl + ^"t^'ruVl + • • • j • (12) 

When n is close to zero the first coefficient on the right hand side is unity. 
We then formally integrate over y^ for a 7^ 1, so that we are left with only one 
integral in the limit n ^ 0: 



exp [~<Ad + Yg't-'.' [i/l - 1) + ■) 

To compute this integral, we can still try to apply a saddle point method 
to the function inside the exponential. This function is proportional to A and 
the first derivative vanishes for two solutions, ya = and yb = 2q(f)rn/ (^0m~3g). 
The values of the second derivatives are equal respectively to minus and plus 
A0^/g. If we expand the integral around ya for example, we obtain 



/m(A) = (^) r dz^zeM-^') + ■■■ (14) 

Vvr V A / 0m JO 

with j;^dz^exp{-z^) = r(3/4)/2 ~ 0.612 708. The density of states is 
finally asymptotically equal to 



2 d 

P(A) = Im-— {exp (-gAlog A + qX - q) 

TT OX 



(15) 



^ Jm (A) exp (2i7rgA 



m] 



m=— 00 
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The approximation (|14|) is not accurate because we have considered only 
one saddle point, but it basically shows that for large m, Im decreases like 
l/yj\m\, so that the series in (|15p is convergent, except for the points where qX 
is an integer. At these points the series diverges, and we obtain peaks in the 
distribution. A more precise computation of (fT^ is to find a path for which 
the integral is convergent. A first transformation is to make the cubic term 
inside the exponential purely imaginary, so that the integral is defined on the 
real axis. This is done by setting z = yi/a such that the coefficient a satisfies 
for example 



«'(-^0^ + |^0D=^ (16) 

The integration path is then along the direction given by the vector a. 
There are 3 solutions for each being proportional to exp(2z7rA;/3), with 
= 0, 1, 2. We will note them ak^m- Now the coefficients of the quadratic term 
in (fT^ is equal to bk^m = ^(Pm^^l mf^l- Therefore Im can expressed as 



with the function 



(17) 



ipa{bk,m) = V2 J dzy/z eiq){—bf^.mz'^) (cosz^ + crsin^;^^ . (18) 

where a = ±1 has to be determined. In order for ipa- to be well defined, 
the real part of bk^m should be positive. For each m, we choose k such as 
Re(6fc^m) > 0. We have checked numerically that there is always only one 
solution k{m) satisfying this condition. Moreover, for a = 1 the series over m 
in ()15|) vanishes numerically, so that a = —1 gives a finite answer. We compute 
numerically ipaibk^m) and /m(A) for every m, and plot in figure 2 the expression 
(fl^j) in the region of localization above A = 1.7, together with data from figure 
1. The asymptotic curve is in agreement for A > 1.9 with the numerical 
results that use direct diagonalization of random matrices over a large number 
of configurations. Divergences occurs each time that q\ is proportional to an 
integer. Near these points the saddle point approximation may not be accurate 
since in the series [equation (|TH)l ] the terms are decreasing slowly with m, and 
we may need more correcting terms in ()13|l . To get a further idea on how to 
simplify the series, we see that, for large A or large m, ak^m behaves like 




exp(2?7rfc/3). (19) 
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This is a good approximation for A ^ exp(y3/4g) 1.21 if g = 20 and 
m = 0. For a ratio between the two terms inside the brackets in (|16|) equal 
to 1/10, with m = 0, we obtain a value A = 1.06. In this approximation, 
bo,m = {9qX/2y^^ is real and positive, and independent of m, the other 
solutions have negative real parts. We then obtain the following behavior 
for large A 

/.„(A) expitarc/A)^ (20) 




We find that in this limit Im is directly proportional to the inverse of (pm- 
The argument of 0^ is 

arg(0^) = ^ arctan ^^^^ + (21) 

with arg(0o) = ''i"/4, and 6{x) is 1 for x > and zero for x < 0. An 
approximation of the probability distribution function is then given by 



exp(— gA log X + qX — q) 




log A I 

sin(27rgAm — arg(0m) + (T7r/4) 



(l + 47r2mV log^ A 



1/4 



f22) 



The term m = gives the monotonic part of the asymptotic curve, and in 
order for this term not to be zero, we have to take a = —1. we can moreover 
replace ipo- by its asymptotic value 

1 2-^'" 




A further approximation is to replace the denominator in (j22p by (1 + 
TT^m^/ log^ A). The resulting function is less accurate than the numerical 
integration of (fT^ . but it should be accurate enough for very large A. Replacing 
the arctan [equation (j2H)] by ±7r/2, the series can be summed up using 
standard formulas. Instead, we would like to study the behavior of the f(t) 
for large times, which is connected to the behavior of P at small argument. 
Indeed, the main contribution of the integral in © comes from very small A, 
and therefore we need the asymptotic behavior of P in this region. We make 
the hypothesis that the saddle point equation [equation should still be 
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valid for both large and small A, since the asymptotic solutions [equation (jSI)] 
in the two cases are the conjugates of each other using the mapping A 1/A. 
We might however modify the formula [equation (|22p] since ip„ has a different 
behavior for small A 



and we also need to define the new arguments of the complex saddle points 



with arg((/)o) = — 7r/4. In that case, a = 1, so that ipi is roughly constant 
for small A [equation (j21I)]- This result show that should still be an 
asymptotic solution for small and non zero values of A, with ipi instead of 
This may explains the oscillations seen in the low eigenvalue region in 
figure 1 (inset), for 0.2 < A < 0.5, with g = 20 and = 800. we can notice 
that, despite the fact that the solutions from (jSJ are conjugate by A — >• 1/A, 
and thus possess some symmetry, the two regions A -C 1 and A 1 appear 
not to have symmetrical asymptotic distributions, since the global distribution 
itself inside the mobility edges is not symmetric. 

A further study would be to precise the physical meaning of these peaks 
in the localized region, and this actually appears in the structure of the real 
eigenvectors (see 0). It may correspond to the particle trap in a region 
of low or high connectivity (defects) as suggested by BM. An interesting 
model given by BM is a Cayley tree with a defect on the central site, with 
a connectivity c instead of g + 1 for the other sites. They found some localized 
states corresponding to, for example, strong c. This localized states can 
disappear if a connectivity c' for the surrounding neighbors is introduced 
and if c' is small enough, giving rise to a screening effect. Also, it is not 
clear however why the structure is 1/g periodic and if the peaks observed are 
Dirac peaks or are diverging with a power law like our result suggests. The 
asymptotic result ()22j) is also a complement to approximations given in the 
BM's work, and also in a different way by reference where the central 
limit theorem is used to computed the fixed-point function of an implicit 
equation giving the probability distribution for an effective Hamiltonian. One 
of the BM's approximations is based on a single defect approximation, which 
consists to allow fluctuations of the connectivity of a single site within an 
effective medium. They were able to find numerically the peaks in the localized 





(25) 
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region and gave the value of the different weights for a given with a good 
approximation. In reference jTHj, peaks seem not to diverge for large A, they 
appear to be rounded, but they are observed with a shift to the right as noted 
by the author. In the small A limit, peaks seem to be sharper, but the central 
limit theorem may not work in this case, even if the position of these peaks are 
correct. Interestingly the correct shape of the distribution in the delocalized 
region is found with high accuracy. 

I thank Remi Monasson for introducing me on this subject and for useful 
discussions. 
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Figure 1. Probability density function P(A) averaged over 4500 samples 
with N=800 and q=20. The inserts show oscillations at the edges of the 
distribution. 




Figure 2. Probability density function P(A) from the asymptotic form 
(|15|l superposed with the numerical diagonalization curve from figure 1 in 
the large A region. 



